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ABSTRACT 

We present a study of the Kelvin-Helmholtz instability in a weakly ionised, multi- 
fluid MHD plasma with parameters matching those of a typical molecular cloud. The 
instability is capable of transforming well-ordered flows into disordered flows. As a re- 
sult, it may be able to convert the energy found in, for example, bowshocks from stellar 
jets into the turbulent energy found in molecular clouds. As these clouds are weakly 
ionised, the ideal magnetohydrodynamic approximation does not a pply at scales of 
aroun d a tenth of a parsec or less. This paper extends the work of I Jones fc Dowries! 
(|2011f ) on the evolution of the Kelvin-Helmholtz instability in the presence of multi- 
fluid magnetohydrodynamic effects. These effects of ambipolar diffusion and the Hall 
effect are here studied together under physical parameters applicable to molecular 
clouds. We restrict our attention to the case of a single shear layer with a transonic, 
but super- Alfenic, velocity jump and the computational domain is chosen to match 
the wavelength of the linearly fastest growing mode of the instability. 

We find that while the introduction of multifluid effects does not affect the lin- 
ear growth rates of the instability, the non-linear behaviour undergoes considerable 
change. The magnetic field is decoupled from the bulk flow as a result of the ambipolar 
diffusion, which leads to a significant difference in the evolution of the field. The Hall 
effect would be expected to lead to a noticeable re-orientation of the magnetic field 
lines perpendicular to the plane. However, the results reveal that the combination with 
ambipolar diffusion leads to a surprisingly effective suppression of this effect. 
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1 INTRODUCTION 

The Kelvin-Helmholtz (KH) instability can occur any- 
where that has a velocity shear and, as a result, is an 
important instability in almost any system involving flu- 
ids. The instability has been studied i n a variety of as- 
trophysical systems, from solar winds jAmerstorfer et al.l 



a ys 

120071 ; iBettarini et al . 2006; Hascgawa et al.ll2004 ) and pul- 
sar winds dBucciantini fc Del Zanna 20061) to thermal flares 
jVenter fc Meintied bOOff T Due to its ability to drive mix- 
and turbulence, the KH instability has been 



mg 



sidered relevant in protoplan etary disks (|Johansen et al.l 



120061 ; iGomez fc Ostrikej 



tospheres ( Li fc Naravanl 



2005) 



2004) 



accretion disks and magne- 
_ , and other jets and outflows 

|Batv fc Keppensll2006ir 

One environment in which the presence of turbulent en- 
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ergy is of particular interest is molecular clouds. One possi- 
ble source of this turbulent energy is the interaction of pro- 
to stellar jets with t he su rrounding cloud, as first proposed 
bv lNorman fc Silkl ([l980). The transfer of momentum from 
a jet to the cloud is most l ikely to occu r through so-called 
"prompt entrainment" (e.g. lDvsonlll984h . As the supersonic 
protostellar jet propagates into the surrounding cloud, it 
forms a bow shock. This shock accelerates molecular cloud 
material, imparting momentum to it. It is worth noting that 
the momentum imparted is relatively well ordered: a further 
process must occur to convert this momentum into turbulent 
motions. 

Early studies of the KH instability in protostellar jets 
were carried out using perturb ative linear analy sis. These 
initial studies (as reviewed by iBirkins haw 1991;) examine 
the fluid equations under set-ups of various combinations 
of magnetic field and shear layers, etc. By approaching 
the problem using a simple mathematical treatment, the 
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effect of each physi cal parameter can be followed closely. 
lHardee et alj (jl997) , for example, investigate the KH in- 
stability in jets by solving the dispersion relations for KH 
modes over a wide range of perturbation frequencies. How- 
ever, these studies are limited to the linear regime of the 
instability, while the evolution of s tellar jets is very much 
governed by nonlinear phenomena (|Bodo et al.| [l994). Sub- 
sequent studies have followed the growth of the KH insta- 
bility in jet s usin g ti me-dependent numeric al simulations. 
IStone et al. | (|l997h and lDownes fc Ravi | 1998), to name just 
a few, impose linear perturbations onto an initially stable 
set-up to observe the behaviour of the instability into the 
non-linear regime and the effect it can have in stellar jets. 

While many authors have investigated the role 
of the KH instability in general m agnetised and un- 



magnetised astrophy s ical flows (e . g . Frank et alj 19961; 
Malagoli et"aH Il99a |Hardee et at] 1 19971 : iDownes fe Ravi 



1998; Kcp pens et al 



1999), these studies have investigated 
the KH instability in the context of either hydrodynamics 
or ideal magnetohydrodynamics (MHD). These assumptions 
are, however, not always valid, particularly in weakly ionised 
systems. For example, in molecular clouds we know that 
non-ideal e ffects are important at length scales below about 
0.2 p c (e.g. lOishi fc Mac Low! l200rj ; bownes fc O'Sullivanl 
12009] ) and hence it is of interest to explore the KH insta- 
bility in the context of either non-ideal MHD or, preferably, 
fully multifluid MHD. 

In recent years the empha sis of KH studi e s has turned 
to including non-ideal effects. iKeppens et alj (1 19991 ) stud- 
ied both the linear growth and subsequent nonlinear satu- 
ration of the KH instability using resistive MHD numerical 
simulations. The inclusion of diffusion allowed for magnetic 
reconnection and non-ideal effects were observed through 
tearing instabilities and the formation of magnetic islands. 
The case in support of using numerical diffusion in order 
to simulate non-ideal MHD e ffects was argued the follow- 
ing year bv lJeong et alj (|2000f ). as analogous to the similar 
practice used to simulate non-ideal hydrodynamic flows of 
high Reynolds number. 

iBirk fc Wiechenl (|2002h examined the case of a partially 
ionised dusty plasma, using a multifluid approach in which 
collisions could be included or ignored. They found that col- 
lisions between the neutral fluid and dust particles could 
lead to the stabilisation of KH modes of particular wave- 
lengths. The unstable modes led to a significant local ampli- 
fication of the magnetic field strength through the formation 
of vortices and current sheets. In the nonlinear regime they 
observed the magnetic flux being redistributed by magnetic 
reconnection. It was suggested that this could be applicable 
to dense molecular clouds and have important implications 
for th e magnetic flux loss problem l|Umebavashi fc Nakanol 
ll990T ). 

A comprehensive study was carried out by IWiechenl 

(2006) which demonstrated the effect of dealing with the 
plasma using a multifluid scheme. This study focused on 
the effect of varying the properties of the dust grains. The 
results of the simulations led to the conclusions that more 
massive dust grains have a stabilising effect on the system 
while higher charged numbers have a destabilising effect. It 
was found that there is no significant dependence on the 
cha rge polarity of the du st. 

IPalotti et al.l (|2008) also carried out a series of simu- 



lations using resistive MHD. They found that, following its 
initial growth, the KH instability decays at a rate that de- 
creases with decreasing plasma resistivity, at least within 
the range of resistivities accessible to their simulations. They 
also found that magnetisation increases the efficiency of mo- 
mentum transport, and that the transport increases with 
decrea sing resistivity. 

In|j ones fc Downes! l|201ll . henceforth Paper I) we ex- 
amine the behaviour of the KH instability in the presence 
of multifluid effects. We found that, while the linear growth 
rates of the instability are unaffected by multifluid effects, 
the non- linear behaviour was remarkably different. The in- 
clusion of ambipolar diffusion leads to the removal of large 
quantities of magnetic energy while the Hall effect, if strong 
enough, proved capable of introducing a dynamo effect. This 
leads to continuing strong growth of the magnetic field well 
into the non-linear regime and a lack of true saturation of 
the instability. 

In this paper we perform a numerical simulation of the 
complete evolution of the KH instability in a weakly ionised, 
multifluid plasma as found in molecular clouds. We address 
the case of a shear layer with a transonic velocity difference. 
We choose a magnetic field strength typical of dense molec- 
ular clouds and this yields a velocity difference across the 
shear layer which is highly super- Alfvenic (see section [2]). 
It is worth noting here that the precise value of the Alfven 
num ber is known to influence the evolut ion of the instability 
(e.g. lJones et al.lll997l : iBatv et a"l]|2003l ) . The work here ca n 
be most directly compared to Case 4 of I Jones et all l|l997h . 
although it should be borne in mind that our boundary con- 
ditions are slightly different. The molecular cloud material is 
simulated by four individual fluids: a neutral fluid, an elec- 
tron fluid, a positively charged metal ion fluid and a fluid of 
large, negatively charged dust grains. The non-ideal effects 
of ambipolar diffusion and the Hall effect are included in the 
simulation, and their effect on the linear development, sat- 
uration and subsequent behaviour of the instability is anal- 
ysed in detail. 

The aim of this work is to investigate the growth and 
saturation of the KH instability under the influence of the 
multifluid effects found in molecular clouds. The KH insta- 
bility is of particular interest as a possible means of convert- 
ing the ordered energy injected into the cloud by protostellar 
jets into the turbulent energy observed. In Paper I we ran 
simulations with parameters chosen to simulate very high, 
medium and very low magnetic Reynolds number systems 
and with parameters chosen to ensure ambipolar-dominated 
flows and Hall-dominated flows in order to develop a full 
understanding of the roles of each. In this paper, we investi- 
gate the combined effects of these two multifluid effects on 
the instability under parameters that describe the specific 
environment of molecular clouds and compare these to roles 
played by each on the development of the instability as stud- 
ied in Paper I. In doing so, we determine the behaviour of 
the instability in a physical application which, in practice, 
should prove observable. 

In section [5] we outline the multifluid equations used by 
the code for this analysis, the physical model being simu- 
lated, and the computational parameters employed. In sec- 
tion 13 we describe how the ability of the code to simulate 
the KH instability has been validated, in both the cases of 
ideal MHD and in the presence of multifluid effects. In sec- 
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tion[4]we detail the results in both the linear and non-linear 
regimes, indicating where the non-ideal behaviour can be at- 
tributed to individual effects of ambipolar diffusion and the 
Hall effect, or the combination of the two. 



2 NUMERICAL SETUP 

The simulations described in this work are perform ed using 
the HYDRA code jO' Sullivan fc DowneslkoOrj . 120071 ) for mul- 
tifluid magnetohydrodynamics in the weakly ionised regime. 
We further assume the flow is isothermal. The assumption of 
weak ionisation allows us to ignore the inertia of the charged 
species and allows us to derive a (relatively) straightforward 
generalised Ohm's law. The resulting system of equations, 
given below, incorporates finite parallel, Hall and Pederson 
conductivity and, as such, is valid in molecular clouds. In 
these regions the viscous lengthscales are much smaller than 
those over which non- ideal effects are important. This leads 
to high Prandtl numbers and plasma flows in these regions 
can be considered to be effectively inviscid. 

Although protostellar jets in molecular clouds tend to 
be hypersonic l|Ravlll988l ). this study would be appropri- 
ate to the KH instability which might arise due to the 
shear between the transonic flow of the swept-up material 
along the edge of the bowshock and the ambient medium, 
or to the shear along the contact discontinuity just behind 
the bowshock. It has been found that the KH instability 
plays a greater r ole in the transon i c, rather than the su- 
personic, regime l|Miural 1 19901 . Il992l ; iKobavashi et al.ll200Sl ; 
iFrank et al.lll996T ). Thus the low Mach number flows, taken 
in concert with the isothermal assumption, means that fea- 
tures in the flow such as shocks are unlikely to create regions 
of high ionisation. 

In this work we compare the results of three simula- 
tions: the first is a full, multifluid MHD simulation of the 
development of the KH instability in a plasma with molec- 
ular cloud properties, the second is of the development of 
the KH instability in the ideal MHD limit, and the third 
is a simulation of the development of the instability in a 
hydrodynamic system. 



JV-l 



OLiPi = 0, 



(7) 



where pi, q^, B, and J are the mass density and velocity of 
fluid i, the magnetic field and current density, respectively. c s 
denotes the sound speed, and Oi and Kio are the charge-to- 
mass ratios and the collision coefficients between the charged 
species i and the neutral fluid, respectively. 

These equations lead to an expression for the electric 
field in the frame of the fluid, E', given by the generalised 
Ohm's Law 

E' = E q +E h + Ea, (8) 

where the components of the field are given by 

Eo = (J-a )a , (9) 

E H = J x a H , (10) 

E A = -(J x a H ) x a H , (11) 

using the definitions ao = /oB, aH = /hB, aA = /aB, 
where fa = y/ro/B, /h = r H /B and f A = ^rX/B. The 
resistivities given here are the Pederson, Hall and ambipolar 
resistivities, respectively, and are defined by 



ro 



va 



9 i 9 ' 

a H + a A 



CA 



where the conductivities are given by 

JV-l 



o"o = — ^ aiiPiPi, 



JV-l 

1 a iP 
°"H = — X 



B f-< l+/3f 

1 — 1 

JV-l 
1 \ ~> OiiPiPi 
°"A = — 



(12) 
(13) 
(14) 

(15) 
(16) 
(17) 



2.1 Multifluid equations 

The code HYDRA solves the following equations for a system 
of TV fluids. The simulations described in this paper consist 
of TV = 4 fluids, indexed by i — for the neutral fluid 
and i = 1,2,3 for the electron, ion, and dust grain fluids 
respectively. The equations to be solved are 



dpi 



+ V-( / o i q i )=0 (0<i^TV-l), 



dt 

dt 

^ + V ■ (q B - Bq ) = -V x E', 



+ V • (pqoqo + c s pl) = J x B, 



(1) 
(2) 
(3) 



ai/9i(E + qiXB) + pipo#to(qo-qi) = (1 ^ i < TV-1),(4) 

V ■ B = 0, (5) 

V x B = J, (6) 



where the Hall parameter Pi for a charged species is given 
by 

ft = (18) 

To solve these equations numerically we use three dif- 
ferent operators: 

(i) solve equations (1), (2), (3), including the restriction 
of equation (5) and for i = 0, using a standard second or- 
der, finite volume shock-capturing scheme. Note that for this 
operator the resistivity terms in equation (3) are not incor- 
porated. Equutjori_^ is_Jncorporated using the method of 
Dedner jDedner et al.ll2002l ). 

(ii) Incorporate the resistive effects in equation (3) us- 
ing super-time-stepping to accelerate the ambipolar diffu- 
sion term and the Hall Diffusion Scheme to deal with the 
Hall term. 

(iii) Solve equations (4) for the charged species velocities 
and use these to update equation (1) with i = 1, . . . , TV — 1. 
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These operators are applied using Strang operator splitting 
in order to maintain the second order accuracy of the overall 
scheme. We refer the reader to lO'Sullivan fc Downed {2006, 
l2007h for a more detailed description. 



2.2 Physical parameters 

This study aims to simulate the multifluid effects on the 
KH instability in the environment of a molecular cloud. In 
a jet-driven bowshock for example, both the swept-up ma- 
terial in the bowshock, as well the surrounding region, are 
understood to be made up of molecular cloud material. For 
this reason, the physical parameters observed in molecular 
clouds are applicable to this system. 



2.2.1 Properties of bulk fluid 

The molecular clo ud material is mode lled after the sin gle- 
size grain model in lWardle fc Nd (|l99gh fsee ljonesll201ll . for 
further details). The plasma is weakly ionised, and therefore 
consists primarily of a neutral particle fluid. The remain- 
der of the plasma is made up of an electron fluid, a posi- 
tively charged metal ion fluid, and a fluid of large, negatively 
charged dust grains. The average mass of the neutral parti- 
cles is taken to be mo = 2.33 m p , where m p is the mass of a 
proton. This corresponds to a fluid of 90% molecular hydro- 
gen and 10% atomic helium by number, which is representa- 
tive of molecular clouds. The ion fluid represents an average 
of ions produced from a number of metal atoms, including 
Na, Mg, Al, Ca, Fe and Ni. These metals have sufficiently 
similar ionisation and recombination rates, and so can be 
modelled collectively as ions of a single positive charge +e 
l|Umebavashi fc Nakanolll9 90). The molecular ions, of which 
HCO + is the most numerous, are significantly less abundant 
than the metal ions, allowing us to neglect them. An average 
mass of m,2 = 24 m p is assigned to the particles of the ion 
fluid, approximately equal to that of a magnesium ion. In 
the single-size grain model, the dust grains are characterised 
by a grain radius r = 0.1 /im and a total grain mass that is 
1% of the total neutral mass. 

The simulations are parametrised by several properties 
that describe the overall system. The sound speed in the 
plasma is calculated from 



7fc B r 

m 



(19) 



where fee = 1.38 x 10 -16 ergK -1 is the Boltzmann constant, 
the adiabatic index 7 is equal to unity for an isothermal sys- 
tem, and the constant temperature is taken t o be T = 30 K, 
which is consistent with molecular clouds (|Wardle fc Nd 
Il999h . This results in a sound speed of 3.26 x 10 cms , 
or c s = 0.326 km s~ . 

The plasma flow in the simulations is required to be 
transonic in order to allow for reasonable growth of the in- 
stability. The sonic Mach number is set as: 

Va 

— = 1, (20) 



so that the relative velocity between the bowshock and am- 
bient medium is equal to the sound speed. In the reference 
frame of the instability, the velocity of the bulk fluid in the 
plasma is given as ±-5- = ±0.188 km s _1 . 



The strength of magnetic fields in molecular clouds can 
vary great ly, but is generally obser ved as a few tens of mi- 
croGauss (fZweibel fc McKeel IT995V The magnitude of the 
magnetic field in this study is chosen as |£?o| = 22.8 pG so 
that the Alfven velocity is: 



VA = 



B 2 

- 0.0326 kms" 1 . 



This gives an Alfven Mach number of 



Ma = — — 

va 



10, 



(21) 



(22) 



so that the plasma is super-Alfvenic. It should be noted 
that there exists a wide range of fluid velocities and mag- 
netic field strengths in molecular clouds, in which the av- 
erage thermal and magnetic energies are approximately in 
equipartition. However, for this study, we have confined our- 
selves to regions of transonic and super-Alfvenic flows as the 
KH instability in a sub-Alfvenic flow wo uld be stabilised by 
the magnetic field l|Chandrasekharlll96ll '). This study is also 
applicable to regions with stronger magnetic fields that are 
not parallel to the plane of plasma flow, but in which the 
strength of the magnetic field as p rojected onto the plane 
of the plasma flow is super-Alfvenic (I Jones et al.ll 19971 ). Our 
choice of a transonic, super -Alfvenic flow is sim ilar, though 
not identical, to Case 4 of I Jones et al.l (ll997T ) and is thus 
not in the "very weak" field regime, but rather the weak 
field regime. In this regime the magnetic field can be ini- 
tially wound up by the growth of the instability and then 
becomes subject to reconnection due to the ensuing field re- 
versals. The final state of the system in this case, at least 
in ideal MHD, is that of a somewhat disordered, widened 
shear layer which is no longer subject to the KH instability 
on wavelengths permitted by the computational set-up. 

2.2.2 Fluid masses and densities 

A typical number de nsity of particles in a molecular cloud is 
n « n = 10 6 cm~ 3 (|Draine et al.lll983l V The number den- 
sity of the remaining fluids c an th en be derived from figure 
2 in lUmebavashi fc Nakanol fl990). The number density of 
the grain fluid is calculated to be n-j, « 1.85 x 10 -14 no. The 
number density of the i on an d electron fluids, according to 
lUmebavashi fc Nakanol (l990), are seen to be approximately 
10 _8 no. However, these parameters are for a molecular cloud 
with a magnetic field of mag nitude 1 mG. I n orde r to ob- 
tain similar conductivities as IWardle fc n3 (|l999l ) within 
the cloud but for a weaker magnetic field, the number den- 
sity of the ion fluid is taken to be n2 ~ 1.02 x 10~ n no. 
The number density of the electron fluid is not set explicitly 
at this point, it is instead defined by the requirement that 
charge neutrality holds at all times. It will be seen that the 
number and mass density of the electron fluid is small, as 
expected. 

The mass densities of the three dominant fluids can then 
be calculated. It is found, for the neutral fluid, that 



po = mono 



3.89 x 10" 18 gcm" 3 . 



For the ion fluid, 

p 2 = m 2 n 2 = 4.08 x 10" 28 gem" 3 . 



(23) 



(24) 



The total mass of the dust grains in molecular clouds is 
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taken to be equiva lent to 1% of the total neutral mass 
jWardle fc NeH l999V This equates to p 3 = 0.01 p = 3.89 x 



10 20 gem 3 . Equivalently, it can be written that 
- = 0.01^1^ = 1.26X10"^. 



(25) 



The charge-to-mass ratios for the three charged fluids 
can be calculated from their charges and masses, as given 
above, and converted into cgs units using lCw 3xl0 9 statC. 
All charged fluids are taken to be singly-charged, and their 
charge-to-mass ratios are thus found to be 



Q2 



a 3 



and 



ax 



±i^ L6X J° °= 1.2 x 1Q» statC g-\ (26) 
mi 24m p 



-le 
m 3 

-le 
mi 



-2.28 x 10 2 statC i 



— = -5.27 x 10 17 statC g _1 . 



(27) 



(28) 



From these values, the necessary mass density for the elec- 
tron fluid can be calculated from 



(29) 

' gem" 3 . This is much smaller than 
the mass density of the other three fluids, as expected. 



axpi + Q2/52 + a-ip-j, = 
to give p 1 = 9.25 x 10" 33 ■ 



2. 2. 3 Properties of charged fluids 

The multifluid effects are implemented in the system 
through collisions between the various fluids. Due to the 
assumption of weak ionisation, only collisions between the 
neutral fluid and each charged fluid are included. Rate co- 
efficients fo r momentum transfe r for each fluid interaction 
are given in lWardle fc Ng| l|l999l ). For collisions between the 
neutral and dust grain fluids, the collision rate is 



-Kdus 



where 



< ov >3 
m3 + mo ' 



(30) 



< av > 3 = nr 2 f^K) 1/2 = 2.i 8 x lO" 8 cmV 1 . (31) 

The collision rate can then be calculated as 
2.18 x 10 _5 cm 3 s _1 



K, 



< OV >3 



3.0 



m 3 + m 1.26 x 10 12 m p + 2.33m p 
= LOxlO^mVV 1 . 



(32) 



Similarly, for the ion and electron fluids, it can be found 
that 

< av > 2 _ 1.6 x 10^ 19 cm 3 s~ 1 
2,0 rri2 + mo 24m p + 2.33m p 

= 3.64 x 10 13 cm 3 g" 1 s" 1 

and finally, following IWardle fc Ng| (|l999l ). 



(33) 



< av >i 



l X 10^ cm 2(^) 1/2 

V y7rm c / 



mi + mo m c + 2.33m p 

= 1.16 X 10 16 cm 3 g _1 B _1 . 



(34) 



Having determined the primary parameters of the flu- 
ids, the secondary parameters can be calculated. The Hall 



parameter is a measure of how well tied a particle is to the 

<*i( — ) TT • ,1 

— — . Using the quan- 

K i,0P0 b M 



magnetic field and is defined by /3; 



tities above, the Hall parameters for each fluid can be cal- 
culated. 



(1.2 x 10 13 )( 



13\ / 2.28X 10" 



(3.64 x 10 13 )(3.89 x 10" 18 ) 



6.43 x 10 1 



P 3 = 



(-2.28xl0 2 )( 2 f x x 1 , 1 ) 'i',7 J ) 
(1.0 x 10 7 )(3.89 x 10" 18 ) 

(-5.27 x 10 17 )( 2 28 X X 1 ;? U " 5 ) 



= -4.3 x 10" 



-8.92 x 10 4 



(35) 



(36) 



(37) 



(1.16 x 10 15 )(3.89 x 10- 18 ) 

The final calculations are to solve for the conductivities. 
From the definitions given in Sect. 12.11 the Pederson, Hall 
and ambipolar resistivities can be calculated as: 

JV-l 

a O = ~^^2 ai PiPi = 413 - 6 S_1 ( 38 ) 



OH 



N-l 

1 \ ^ an Pi 



°~A = — r 



CtipiPi 



b ^ i + m 

■ i 1 



0.01 s" 1 



0.1s" 



(39) 



(40) 



It can be seen that thes e values are equal to those calculated 
bv lWardle fc NtJ (|l999l ) for a weakly ionised molecular cloud. 

It should be noted that this system consists of three 
charged fluids, while our multifluid KH study in Paper I con- 
sisted of only two. In the previous study, the ion fluid was 
responsible for both the ambipolar and the Hall conductiv- 
ity. With three charged fluids in this study, the ion fluid is 
now responsible only for the ambipolar conductivity, while 
the dust grain fluid provides the largest contribution to the 
Hall conductivity. This results in a negative value for the the 
Hall conductivity, as the dust grains are negatively charged. 
The only consequence of this is that the re-orientation ex- 
perienced by the magnetic field as a result of the Hall effect 
is now in the opposite ^-direction to before. 



2.3 Normalised computational parameters 

The physical parameters for molecular clouds, outlined 
above, are transformed into dimensionless units before be- 
ing passed to the numerical code. For this purpose, practical 
characteristic length, time and mass scales must be deter- 
mined. 

The lengthscale of the system is first chosen. The pro- 
tostellar jet is understoo d to be approximately a parsec or 
two in length (|Ravlll987f ). It is chosen to focus on a fraction 
of this length, opting to simulate a region along the edge of 
the bowshock of length 0.2 pc. It has also been found that 
multifluid effects play an im portant role at this lengthscale 
l|Downes fc O'Sullivajj|200gh . The lengthscale is thus set to 



[L] = 0.2 pc = 6.17 x 10 cm. 



(41) 



The timescale is chosen such that the sound speed is of order 
unity. The sound speed has been calculated as c s = 3.26 x 
10 4 cms _1 . This implies a time scale 
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[T] =1.9 x 10 13 : 



(42) 



Using these units, the velocity of the plasma flow is now 

(43) 



,V 
± — v 
2 y 



± 1 -y. 

2 y 



Finally, a mass scale is chosen. This is set so that the 
mass density of the plasma is of order unity. The mass den- 
sity has been calculated as p — 3.89 x 10~ 18 gcm -3 . Us- 
ing the length scale chosen above, this gives a characteristic 
mass scale 



[M] = 9.15 x 10 3 



(44) 



The magnitude of the magnetic field, Bo = 22.8 fiG, can 
also be transformed into a dimensionless quantity, by 

mm 1 ' 2 



IB 



Bo 



0.1 



(45) 



v/4^ [M]Va 

From this calculation onwards, the units of B, E and 
J are written such that the factors of 4n and the speed of 
light no longer appear. 

It is possible now to transform the conductivities, and 
the resulting resistivities, into dimensionless units. The fol- 
lowing values are found: 

ro = — = 6.2 x 10~ 9 



th = 



on 



= -0.00352 



0.0352 



(46) 
(47) 



(48) 



These values can be seen to correspond to those imple- 
mented in Paper I. The amount of ambipolar resistivity 
is equivalent to the simulation with high ambipolar resis- 
tivity in Paper I (ambi-high-hr) with magnetic Reynolds 
number Re m = 28.4, while the amount of Hall resistivity 
is equivalent to the simulation with moderate Hall resistiv- 
ity in Paper I (hall-med-hr) with magnetic Reynolds number 
Re m = 284. 



2.4 Grid initialisation and properties 

As in Paper I, these simulations are carried out on a 2.5 D 
slab grid in the xy-plane. The initial set-up used was that 
of two plasmas flowing anti-parallel side-by-side on a grid of 
size x — [0, 32L] and y — [0, L]. The plasma velocities are 
given by +-£■ and — ^ in the y-direction. The unit of time 
used is the sound crossing time, i s = — , the time taken for 
a signal travelling at the speed of sound to cross the grid in 
the y-direction. 

We use periodic boundary conditions at the high and 
low y boundaries. Since we wish to study not only the initial 
growth phase of the instability, but also its subsequent non- 
linear behaviour we must ensure that waves interacting with 
the high and low x boundaries do not reflect back into the 
domain to influence the dynamics. Previous test simulations 
for various parameters have shown that a large width of 32L 
is necessary to ensure this. We use gradient zero boundary 
conditions at the high and low x and z boundaries. 

The grid is chosen therefore to consist of 6400 x 200 x 1 
cells, in the x, y, and z directions respectively. This resolu- 
tion was chosen on the basis that it reproduces the initial 



linear growth of the ideal MHD system in iKeppens et al.l 
(1999). Resolution studies were performed to confirm the 
resolution as being appropriate for multifluid MHD as well 
(see Sect.GO). 

The plasma velocity field is initiated with a tangential 
shear layer of width 2a at the interface at x — 16L. This 
velocity profile is described by 



v = — tanh I - 



16L\ . 



(49) 



The width of the shear layer is chosen to be j- — 0.05, or 
approximately 20 grid zones. The magnetic field is initially 
set to be uniform and aligned with the plasma flow. 

The initial background for all four fluids in the system 
is now an exact equilibrium. The initial velocity field, Vo is 
then augmented with a perturbation given by 



Sv x = 5Vq sin(—k y y) exp 



(x - 16L) S 



(50) 



where SVo is set to 10 -4 Vo. The wavelength of the per- 
turbation is set equal to the characteristic length scale, 
A = jp = L, so that a single wavelength fits exactly into 
the computational domain. This maximises the possibility 
of resolving structures t hat are small relat ive to the initial 
perturbed wavelength, (|Frank et al.l Il996l ). The perturba- 
tion attenuation scale is chosen so that it is larger than 
the shear layer, but small enough so that the instability 
can be assume d to interact only minimally with the x- 
bou ndaries (see Palotti et al.ll2008 | ). and is set using J = 0.2 
fsee IKeppens et al.l Il999l : iPalotti et alj l2008h . Finally, the 



wavenumber k y is chosen to b e 2n in order to maxi mise the 
growth rate of the instability (jKeppens et al.|[l999l ). 

We note that our choice of periodic boundary conditions 
on the y boundaries limits the behaviour of the system to 
some extent: wavelengths longer than the grid domain can 
not gro w and coalescence of large-scale vortices will not oc- 
cur (e.g. iBatv et al-lkoOSl ) . Furthermore, our choice of gradi- 
ent zero boundary conditions in the x direction implies that 
we are studying only the surface modes of the instability. 
This latter point, however, is not a major restriction as the 
body modes will not be important in, for example, the de- 
velopment of the KH instability in the bowshocks created 
by young stellar outflows. 



3 VALIDATION OF NUMERICAL APPROACH 

3.1 Validation of instability growth in ideal MHD 

The set-up described above has been shown in Paper I to 
be valid for producing the KH instability under the ideal 
MHD approximation through comparisons with previously 
published literature. Both hydrodynamic and ideal MHD 
simulations were run using the HYDRA code. Comparisons 
between the linear growth of the simulated hydrodynamic 
instability and the growt h rate calculated analyti cally for 
the same wavenumber bv lMiura fc Pritchettl (| 19821 ) showed 
exceptional agreement. In the non-linear regime, the max- 
imum reached by the mag netic energy in the i deal MHD 
simulation matches that of iMalagoli et al.l (|l996) to within 
10%. This allows us to be confident of the behaviour of HY- 
DRA in simulating the KH instability. 
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3.2 Validation of instability growth in multifluid 
MHD 

The inclusion of multifluid effects introduces new length 
scales into the system. These include the diffusion length 
scales of the magnetic field due to ambipolar resistivity, 
and the rather computationally challenging whistler waves 
arising from the Hall effect. The Hall term is handled by 
HYDRA using the explicit Hall Diffusion Scheme (HDS) 
lO'Sullivan fc DowneslfeoOrj , 120071 ). Although the code nat- 
urally does not resolve waves of vanishing wavelength, reso- 
lution studies were performed in Paper I of both ambipolar 
and Hall-dominated flows and the results indicated that a 
resolution of 6400 x 200 x 1 is sufficient to capture the initial 
growth and saturation of the instability. Subsequently, the 
dynamics are captured at least qualitatively. 

In order to ensure that the smallest-scale dispersive ef- 
fects were in place when examining whether they are suf- 
ficiently resolved, the highest values for the ambipolar and 
Hall resistivity used in Paper I were implemented for these 
resolution studies. As indicated in Sect. 12.31 this highest 
value of ambipolar resistivity in Paper I is equivalent to the 
value implemented in this paper (with magnetic Reynolds 
number i?e m = 28.4), while the amount of Hall resistivity 
implemented in this study is equivalent to the simulation 
with only moderate Hall resistivity from Paper I (with mag- 
netic Reynolds number i?e m = 284). 

We are, therefore, confident that the multifluid dynam- 
ics resulting from the non-ideal effects in these simulations 
are well resolved and that our conclusions as to the physical 
processes occurring are well-founded. 



4 RESULTS 

In order to observe the differences in the evolution of the 
KH instability with the inclusion of multifluid effects, sev- 
eral aspects of the evolution are examined and compared 
to those from simulations carried out in ideal MHD and in 
pure hydrodynamics (see |Jonesll201ll . for further details of 
the simulation results). 

The study of the growth of the instability is carried out 
through measuring the evolution of a number of parameters 
with time. In particular, we measure the transverse kinetic 
energy 



Eh 



^pvl dx dy 



and the magnetic energy 

E b = i J J { [Bl + B 2 y + Bl] - B 2 } dxdy 

in the system where -Bo is the magnitude of the magnetic 
field at t = 0. Any growth of E^ x is due to the growth of 
the instability, as the entire plasma flow is initially in the 
y-direction, with only a very small perturbation in the x- 
direction, 



4.1 Linear regime 

The KH instability leads to an interaction between the two 
plasmas on either side of the initial interface. In particu- 
lar, in ideal MHD, the plasmas are seen to wind-up, leading 
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Figure 1. Plot of the magnitude and vector field of the bulk 
velocity field for the multifluid MHD case, at time of saturation 
of the instability, t = 8t B . 




Figure 2. The log of the transverse kinetic energy in the sys- 
tem is plotted against time for the ideal MHD (thin line) and 
multifluid MHD (thick line) cases. The dashed line represents the 
hydrodynamic case. There is very little difference between the 
linear growths for the ideal and multifluid cases. 



to the "Kelvin's cat's eye" vortex. The inclusion of multi- 
fluid effects can affect this evolution in a number of ways. 
However, it can be seen that this multifluid set-up does not 
prevent the development of the classic vortex in the neutral 
fluid (see Fig. [[} . 

As described in Paper I, analysis of the transverse ki- 
netic energy in the system allows for study of the growth 
rate of the instability. Figure [2] plots the growth of the ki- 
netic energy resulting from the instability in the ideal and 
multifluid MHD cases, as well as the hydrodynamic case. It 
can be seen that the linear growth of the instability differs 
very little with the inclusion of multifluid effects. It will be 
seen that the instability is developing in a way very different 
to the ideal MHD case, yet the linear growth rates are found 
to be within 1% of each other. 



4.2 Non-linear regime 

4-2.1 Diffusion 

In the ideal MHD case, the plasma is wound up by the KH 
instability and the magnetic field experiences a similar wind- 
ing force as a result of the frozen-in approximation (see Fig. 
[3] upper panel). 

It can clearly be seen that the magnetic field under- 
goes a very different evolution when multifluid effects are 
included (see Fig. [3] lower panel). The inclusion of ambipo- 
lar resistivity into the system allows for decoupling between 
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Figure 3. Plot of the magnitude and vector field of the magnetic 
field for the ideal MHD (upper panel) and multifluid MHD cases 
(lower panel) , at time = 8 t B . In the multifluid MHD case the 
magnetic field no longer demonstrates a wind-up as seen in the 
ideal MHD case, as a result of decoupling and diffusion. 
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Figure 4. The perturbed magnetic energy in the system is plot- 
ted against time in the ideal MHD case (thin line) and multifluid 
MHD case (thick line). The introduction of multifluid effects sig- 
nificantly dampens the amplification of the magnetic energy. 



the various fluids. This breaks the frozen-in approximation 
of ideal MHD. As a result, the magnetic field is able to dif- 
fuse with respect to the bulk fluid. This ambipolar diffusion 
is the source of the altered magnetic field configuration ob- 
served. 

The changes in the magnetic field development can be 
analysed in a more quantitative manner using the plot in 
figure [4] The thin line shows the amplification experienced 
by the magnetic field through the wind-up it undergoes in 
the ideal MHD case. The thick line shows that this amplifi- 
cation is significantly reduced in the presence of ambipolar 
diffusion. Similar results were observed in Paper I in the 
ambipolar-dominated simulation. As this was not observed 
in the Hall-dominated simulations in Paper I, we can deduce 
that this is solely as a result of the ambipolar resistivity. 



4-2.2 Charged fluids 

We know that the introduction of ambipolar resistivity has 
allowed for decoupling of the magnetic field from the neutral 
fluid. We now examine the behaviour of the charged fluids 
themselves. The exact behaviour of each charged fluid can 
be understood by examining its density profile and velocity 
field during the development of the instability. The state of 
each of the four fluids in the system has been plotted in 
figure [S] at the time of saturation of the instability. 

It can be seen that the mass density of the dust grain 
fluid closely reflects that of the bulk fluid, signifying a strong 
coupling between the two. This would be expected due to its 
relatively low Hall parameter (see equation l36p . On the other 
hand, the ion and electron fluids more closely reflect the 
configuration of the magnetic field, implying that they are 
still strongly coupled to the magnetic field lines, as expected 
by their high Hall parameter (equations I35I and I37[) . The 
decoupling of the ion and electron fluids from the neutral 
fluid is the source of the ambipolar diffusion in the system. 
The magnetic field is tied to the neutral fluid only through 
the coupling of the charged fluids with the neutrals, so a 
low collisional coupling between the charged fluids and the 
neutrals allows for the magnetic field to diffuse relative to 
the bulk fluid. 

The various dynamics discussed above are confirmed in 
figure[U] These plots of the transverse kinetic energy for each 
of the four fluids clearly demonstrate the behaviour of each. 
We can see that the bulk fluid undergoes further wind-up 
in the multifluid MHD case than in the ideal MHD case 
(see the top panel of figure [6}. This is due to two distinct 
phenomena (see Paper I for more details). In general, the 
system is prevented from as strong a wind-up as seen in 
the hydrodynamic case by the presence of a magnetic field. 
In multifluid MHD, there are two effects at work that limit 
the effectiveness of the magnetic field in suppressing this 
wind-up. Firstly, the introduction of even a small amount 
of ambipolar diffusion causes the magnetic field to experi- 
ence a significant reduction in its amplification. The result- 
ing weaker magnetic field allows the bulk fluid to undergo a 
stronger wind-up. Secondly, with higher levels of ambipolar 
diffusion being introduced into the system, the bulk fluid 
becomes further decoupled from the magnetic field, further 
reducing its effectiveness in opposing the wind-up. 

On the other hand both the electron and ion fluids ex- 
perience a decoupling from the neutral fluid. As they are 
still well tied to the magnetic field, and the magnetic field 
no longer winds up in a manner similar to the bulk fluid, the 
charged fluids undergo less wind-up than in the ideal MHD 
case. This is demonstrated by the plots of the transverse 
kinetic energy of these fluids (figure [5] second and third 
panels), and the lower maxima reached. Finally, the plot 
representing the dust grain fluid (figure [6] bottom panel) 
shows similar behaviour to that observed in the bulk fluid. 
This confirms that the dust grains remain well-coupled to 
the neutral fluid, and are less affected by the magnetic field 
than the other charged fluids. 



4-2.3 Combination of multifluid effects 

The results detailed above can be attributed to the introduc- 
tion of ambipolar diffusion into the system. However, there 
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Figure 5. Plot of the distribution of the mass density, as well 
as the velocity vector field, of the bulk fluid (top panel), electron 
fluid (second panel), ion fluid (third panel) and the dust grain 
fluid (bottom panel) in the multifluid case at time = 8t s . 



is a significant amount of Hall resistivity included in the set- 
up as well. While the system remains ambipolar-dominated, 
results from Paper I indicate that there should be observable 
consequences of including this moderate amount of Hall re- 
sistivity. Most notably, a small but significant re-orientation 
of the magnetic field is expected, from the xi/-plane into the 
z-direction. Plotting the evolution of the magnetic energies 
in each of the three directions however, shows that there is, 
in fact, no significant growth in the z-direction (see figure 
[7J. This is in distinct opposition with the results from Pa- 



Figure 6. The transverse kinetic energy is plotted against time 
for the ideal MHD (thin line), multifluid MHD (thick line) and hy- 
drodynamic cases (dashed line), for the entire system (top panel), 
electron fluid (second panel), ion fluid (third panel) and dust grain 
fluid (bottom panel). 



per I, in which significant growth of the magnetic energy 
in the z-direction was observed for the same level of Hall 
resistivity. 

It is important to recall that the strength of the Hall 
effect depends on the current in the system. The current, in 
turn, depends on the charge densities of the three charged 
fluids. The charge density of the dust grain fluid is a^pz « 
-8 x 10" 18 state cm . However, for the ion and electron 
fluids, the charge densities are much higher, approximately 
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Figure 7. The total perturbed magnetic energy (solid line) in 
the system is plotted against time in the multifluid MHD case. 
This is almost identical to the magnetic energy in the xj/-planc 
(dashed line), while there is very little growth in the magnetic 
energy in the z-direction (dot-dashed line). 



±5 x 10~ 15 statC cm -3 . Therefore the current in the system, 
J = y^ . OiPiVi, is primarily due to the velocity difference 
between the ion and electron fluids. 

Therefore, while the decoupling of the dust grain fluid 
from the magnetic field provides a high Hall resistivity, the 
strength of the Hall effect is in fact also critically depen- 
dent on the dynamics of the ion fluid relative to the electron 
fluid. As the ion fluid Hall parameter is much greater than 1 
it will be somewhat decoupled from the neutral fluid, instead 
being more strongly tied to the magnetic field. This means 
that the relative velocity between the ion and electron fluid 
will be rather small. As a result, the majority of current in 
the system remains parallel to the magnetic field, and the 
strength of the Hall effect, which is proportional to J x B, 
is very small. This accounts for the growth of the magnetic 
field in the z-direction being much less than would naively 
be expected. Put simply, the introduction of ambipolar re- 
sistivity changes which fluids are coupled to which and, in 
doing so, inhibits the Hall effect. 

On a side note, this can also occur in a system with 
only two charged fluids. In this case, it is typically the ion 
fluid that is the larger contributor to both the Hall and 
ambipolar conductivity. The decoupling of the ion fluid from 
the magnetic field leads to Hall resistivity, as any velocity 
difference between the ion and electron fluid can give rise to 
a current with a component perpendicular to the magnetic 
field. However, the ambipolar resistivity arises due to a low 
collision rate between the ion fluid and the neutral fluid. As 
a result, the ions do not, in fact, behave in the same way as 
the neutrals. Their coupling with the magnetic field, though 
weak, is sufficient to ensure only a minimal relative velocity 
between the ions and electrons. In this way, the impact of 
the Hall effect is similarly reduced. 

In a Hall-dominated flow, we would expect to see not 
only a growth in magnetic field strength in the z-direction, 
but also a resulting growth of kinetic energy in the same 
direction, as the plasma is influenced to travel out of the xy- 
plane. With even moderate Hall resistivity, this kinetic en- 
ergy in the z-direction can become comparable to that in the 
^-direction (see Paper I). However, in this multifluid case, 
the magnetic field experiences only a minimal re-orientation 
into the z-direction, and therefore only the fluids that are 
tightly coupled to the magnetic field will experience any no- 
ticeable dynamics in this direction. Figure [8] shows that the 
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Figure 8. The growth of the kinetic energy of the entire system 
(upper panel) and the ion fluid (lower panel) is plotted against 
time in both the ^-direction (solid line) and z-direction (dashed 
line). The growth in the z-direction for the bulk flow is clearly 
negligible, while the growth in the z-direction for the ion fluid is 
very small, but non-negligible. 



bulk flow demonstrates negligible growth of kinetic energy in 
the z-direction (upper panel), while the ion fluid, being more 
closely tied to the magnetic field, demonstrates a very small, 
but non- negligible, growth in this direction (lower panel). 



4.3 Subsequent behaviour 

The KH instability is seen to undergo a very different evo- 
lution in the presence of multifluid effects. In the ideal case, 
the initial wind-up of the velocity field has the effect of also 
winding up the magnetic field, due to strong coupling be- 
tween the two. As the KH vortex is caused to stretch and 
expand by the amplified magnetic field, it reaches the peri- 
odic {/-boundaries of the simulated grid. As the neighbouring 
vortices merge, numerical viscosity allows for magnetic re- 
connection, which results in the creation of magnetic islands 
and secondary vortices. The generation and decay of these 
vortices results in further periods of growth in the trans- 
verse kinetic energy in the system corresponding to periods 
of decreasing magnetic energy and vice- versa. This is clearly 
demonstrated in figure [9] where the peaks in each energy 
are seen to correspond to the troughs in the other. This be- 
haviour continues in the ideal MHD system until the system 
has reached a somewhat disordered state (see figure [10}. 

In the multifluid case, on the other hand, the subsequent 
evolution of the KH instability is quite different. The high 
ambipolar resistivity causes the magnetic field to experience 
very little wind-up, through decoupling and diffusion. As a 
result, the KH vortex undergoes little or no stretching or 
expansion as a consequence of the magnetic field. Following 
the saturation of the instability, the magnetic field is seen 
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Figure 9. The growth of the perturbed magnetic energy (dashed 
line) is plotted against time in the ideal case. Also plotted is the 
transverse kinetic energy (solid line) of the ideal case. It can be 
seen that, after the initial saturation point, the growth of mag- 
netic energy is linked to the decay of kinetic energy, and vice- 
versa. 



Figure 11. The evolution of the transverse kinetic energy (solid 
lines) and magnetic energy (dashed lines) for the ideal MHD case 
(thin lines) and multifluid case (thick lines). It can be seen that, 
after the initial saturation point, the magnetic field has experi- 
enced little growth, and is unable to inject energy back into the 
system. 
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Figure 10. Plot of the magnitude and vector field of the velocity 
field for the ideal MHD case, at time = 20 t B . The system has 
become somewhat disordered. 



to steadily return to its original configuration, and neither 
the velocity field nor the magnetic field undergo a second 
period of growth (see figure ITT]) . Instead, the magnetic field 
eventually stabilises the velocity field, and the plasma re- 
turns almost to its original state, with simple laminar flow 
on either side of a much wider shear layer (see figure [12)) . 
This wider shear layer is no longer conducive to the growth 
of the KH instability. 

As the vortex is broken down and the system returns to 
a stable state, the energy consumed by the instability flat- 
tens out. The system has lost some energy during the initial 
growth of the instability due to ambipolar diffusion remov- 
ing some magnetic energy, but this levels off at later times. 
This is a deviation from the ideal MHD case, in which the 
system continues to lose energy as the KH vortex is broken 
down through reconnection into disordered decay. Figure [131 
plots the total energies of the ideal and multifluid MHD 
cases, as well as the hydrodynamic case. In the hydrody- 
namic case, the KH vortex remains indefinitely, as there is 
no magnetic field to lead to its decay. For these reasons, the 
total energies in the multifluid and hydrodynamic cases are 
seen to level off, while the ideal MHD system continues to 
lose energy through decay. 
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Figure 12. Plot of the magnitude and vector field of the velocity 
field for the high ambipolar case, at time = 20 t B . The system has 
returned to almost laminar flow. 



5 CONCLUSIONS 

A study of the KH instability in a molecular cloud plasma is 
carried out. We restrict our attention to the case of a single, 
transonic, super- Alfvenic shear layer and, as such, the body 
modes of the KH instability are not examined. The KH in- 
stability, in particular the growth of its surface modes, is 
severely reduced for high Mach number flows, while sub- 



0.02 




Figure 13. The total energy in the system of both the ideal 
MHD (thin line) and multifluid MHD (thick line) systems are 
plotted against time. While the ideal case continues to lose energy 
throughout the simulation through disordered decay, the levelling 
off of the energy in the multifluid case indicates that it has reached 
a stable state. 
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Alfvenic flows tend to be stabilised by the presence of the 
magnetic field. Our choice of an Alfven number of 10 for 
the system is realistic in terms of that present in a molec- 
ular cloud and means that the growth of the instability, in 
the ideal MHD approximation where there is nearly no re- 
sistivity, will roll up the magnetic field eventually leading 
to a situation in which tearing mode inst abilities can reor- 
ganise the topol ogy of the magnetic field (I Jones et al.ll 19971 ; 
iBatv et al.ll2003h . We utilise periodic boundary conditions in 
the longitudinal direction and, while a very usual approach, 
it should be remembered that the so-called "inverse cascade" 
observed bv lBatv et all (|2003h will not occur in this system. 

The combined effects of ambipolar diffusion and the 
Hall effect present in such a plasma are examined with ref- 
erence to an ideal MHD system. As molecular clouds are 
dominated by ambipolar diffusion, it is expected that the 
results will closely reflect those of the ambipolar-dominated 
simulation in Paper I and this is indeed observed. However, 
one's naive expectation would be that the Hall effect should 
also be observed through twisting of the magnetic field lines. 
In the initial linear regime, the growth rates of the instabil- 
ity remain unchanged by multifluid effects. As the instability 
develops, ambipolar diffusion leads to less amplification of 
the magnetic energy in the system and therefore a stronger 
wind-up of the neutral fluid. This is very similar to the re- 
sults observed in the ambipolar-only system studied in Pa- 
per I. Subsequent behaviour of the instability, including the 
breakdown of the KH vortex and its return to a stable state 
and laminar flow, also closely reflects that observed in the 
ambipolar-only system. 

A noteworthy result of this study is the lack of the ex- 
pected impact on the system due to the inclusion of the Hall 
effect. This appears to be due to the presence of ambipolar 
diffusion: i.e. ambipolar diffusion inhibits the impact of the 
Hall effect. The set-up presented here includes three charged 
fluids (ions, electrons and dust grains) as well as the neu- 
tral fluid. However, this same result can be observed in a 
system with two charged fluids (ions and electrons) as well 
as the neutral fluid. In the case of two charged fluids, the 
Hall parameter of the ions will be much lower than that of 
the electrons, thereby introducing Hall resistivity. Ambipo- 
lar resistivity then arises when the ion fluid also has a low 
collision rate with, or more properly momentum transfer to, 
the neutral fluid. As the ion fluid is subjected to fewer col- 
lisions with the neutrals, it is the magnetic field that plays 
the dominant role in determining its dynamics, despite the 
relatively weak coupling. This minimises the impact of the 
Hall effect in the system. In the case of three charged flu- 
ids, Hall resistivity arises in the system from the dust grain 
fluid being weakly coupled to the magnetic field. Ambipolar 
diffusion will arise when the collision rate of the ion fluid 
with the neutrals is low. However, when the dust grain fluid 
has a much smaller charge density than that of the ion and 
electron fluids, the effective current is due mainly to the rel- 
ative velocities of the ion and electron fluids and not the 
dust grain fluid. As both the electron and ion fluids remain 
tightly coupled to the magnetic field, their relative veloc- 
ity will continue to be in the direction of the magnetic field 
lines, resulting again in a suppression of the Hall effect in 
the system. 

We conclude that for the parameters observed in 
ambipolar-dominated molecular clouds, the evolution of the 



KH instability is indeed dominated by the presence of am- 
bipolar diffusion. However, the system does not demonstrate 
the expected twisting of the magnetic field lines which would 
arise from the inclusion of the Hall effect as naively expected. 
This suppression of the Hall effect by ambipolar diffusion 
is given an original and detailed explanation based on ba- 
sic physics principles and for plasmas involving a variety of 
charged fluids. 

Finally, we note that it would be fascinating, though 
computationally very expensive, to perform very large scale 
simula tions in which the inverse cascade seen bv lBatv et al.l 
can transfer energy from length scales on which mul- 
tifluid effects are important right up to length scales on 
which ideal MHD is an appropriate approximation in order 
to study the nature of this cascade in detail. For the work 
presented here, however, our chosen length scales of ~ 0.2 pc 
are appropriate for the study of the effect of this instability 
on the bowshocks caused by young stellar outflows. 
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